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We investigate the electronic and other properties of the hypothetical compound LiFeSb in relation 
to superconducting LiFeAs and FeSe using density functional calculations. The results show that 
LiFeSb in the LiFeAs structure would be dynamically stable in the sense of having no unstable 
phonon modes, and would have very similar electronic and magnetic properties to the layered Fe 
based superconductors. Importantly, a very similar structure for the Fermi surface and a spin density 
wave related to but stronger than that in the corresponding As compound is found. These results 
are indicative of possible superconductivity analogous to the Fe-As based compounds if the spin 
density wave can be suppressed by doping or other means. Prospects for synthesizing this material 
in pure form or in solid solution with FeTe are discussed. 
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The finding of high temperature superconductivity 
[Tc ~26K) in electron-doped LaFeAsOi^xFa:,- has re- 
sulted in widespread interest and exploration of related 
materials, some of which have Tc exceeding 55K. In 
particular, superconductivity has been found in iron 
based oxy-arsenides by replacing La with other rare- 
earth metals ^i^i^i^i^ as well as oxygen-free arsenides such 
as doped BaFezAsar^i^ SrFeaAsa^ CaFeaAsa^iii and 
LiPeAsi^ii^iii^. The common structural feature of this 
family of materials is the appearance of Fe-As layers. 
These consist of an Fe square planar sheet tetrahedrally 
coordinated by As atoms from above and below. In addi- 
tion, superconductivity occurs in doped LaFePO j^^i^^i^^ 
although with a lower Tc and in PbO structure a- 
FeSei-a i^^i^^'^" These latter compounds also feature an 
Fe square lattice and a tetrahedral coordination of the 
Fe, though not with As. Importantly, the critical tem- 
perature of FeSei_a; increases strongly with either Te 
substitutioi*i^ or pressure, reaching 27 K,— This high 
value of Tc under pressure implies a relationship with 
the Fe-As superconductors, which is also supported by 
similarities of the properties and theoretical studiesj^ 
At present there is strong interest in finding new high 
temperature Fe-based superconductors and especially in 
finding materials with higher critical temperature. 

One obvious direction is to examine antimonides. This 
is motivated by the fact that the properties of LaFePO 
and LaFeAsO appear to be closely related, suggesting a 
similar mechanism of superconductivity, and furthermore 
the compound with the heavier pnictogen (As) has the 
higher when doped. However, this is highly non-trivial 
from a chemical perspective because Sb has a strong ten- 
dency to form Sb-Sb bonds in compounds. This leads to 
a strong tendency for transition metal compounds of Sb 
to contain more Sb than transition metal, as for exam- 
ple in skutterudite CoSbs and LaFe4Sbi2 or marcasite 
structure FeSb2, although FeSb is a known phase. One 
way forward is provided by noting the structural similar- 
ity of LiFeAs with PbO structure FeSci-^: and FeTci^^:. 



The chalcogenides, whose chemical formulas should more 
correctly be written as Fei+jjSe and Fci+^-Te, occur in a 
tetragonal structure with spacegroup PA/nmm similar 
to LiFeAs, and consist of an Fe square lattice tetrahe- 
drally coordinated with Se/Te ions, the same as in the 
structure of the Fe-As superconductorsi^^i^iiS^i^ These 
chalcogenides form with excess Fe, which occurs in a par- 
tially filled 2c site, in particular the cation site forming 
formally an enlarged tetrahedron around the Fe and ap- 
proximately five-fold coordinated by Te. ^'^i^^'^^ This is 
the same site that is occupied by Li in LiFeAs. There- 
fore, there is a close structural similarity between LiFeAs 
and the chalcogenides Fci+^jSe and Fei+^Te. In partic- 
ular the structure of LiFeAs is obtained by allowing full 
filling of the 2c cation site with Li+ and replacement 
of Te^~ by As^~. Therefore we focus on hypothetical 
LiFeSb since it may be possible to form it, or at the very 
least some range of solid solution between Fei+^^Te and 
LiFeSb should be experimentally accessible, especially 
considering that alloys of related phases containing Te 
and Sb typically form as in e.g. the Bi-Sb-Te, AgSbTe- 
PbTe and AgSbTe-GeTe thermoelectrics, and also that 
there are many known Zintl type phases based on Li, Sb 
and metal atoms. 

The crystal structure of LiFeSb is assumed to be 
isostructural with LiFcAs with the space group of 
P4/nTOmJ^ii^ii^ As shown in Fig. (TJ the Fe-Sb layers 
formed by edge-shared tetrahedral FeSb4 units are alter- 
nately spaced along the c-axis direction, and intercalated 
with Li. Note that the Li, which occurs as Li+, is coordi- 
nated by Sb. The structural parameters were calculated 
by local density approximation (LDA) total energy min- 
imization with the full-potential linearized augmented 
plane wave (LAPW) method.2'^ The calculated tetrag- 
onal lattice parameters are a — 4.0351 A, c = 6.3712 A, 
and internal coordinates Li(2c) (0.25,0.25,0.697), Fe(2a) 
(0.75,0.25,0), Sb(2c) (0.25,0.25,0.228). The Fe-Sb bond 
length is 2.486 A, slightly larger than 2.4204 A for 
LiFeAs/"* which might be attributed to the larger size 



of the Sb^~ anion relative to As"^". The Fe-Fe distance 
is 2.853A, also a bit larger than the corresponding value 
(2.68O9A) in LiFeAs, but still short enough for direct Fe- 
Fe interaction. 

The electronic structure and magnetic properties cal- 
culations were performed within LDA-LAPW method, 
similar to previous reports.— i^i^^ LAPW sphere radii of 
1.8 oq, 2.0 oq, and 2.1 oq were used for Li, Fe and Sb, 
respectively. Converged basis sets were used. These con- 
sisted of LAPW functions with a planewave cutoff deter- 
mined as Rhikmax—&-0 plus local orbitals both to relax 
linearization and to include the semicore states. The zone 
sampling for the self-consistent calculations was done us- 
ing the special k-points method, with a 16x16x8 grid. 
Finer grids were used for the density of states and Fermi 
surface. The lattice dynamical properties were calcu- 
lated through the frozen phonon method^ (or small dis- 
placement method'^^). The required forces were obtained 
through the projector augmented-wave (PAW) method'^^ 
in VASP code, within the generalized gradient approx- 
imation of Perdew, Burke and Ernzerhof (PBE);2^ We 
also fully relaxed crystal structure and calculated elec- 
tronic structures with PEE- PAW method and the results 
show excellent agreement with those by LDA-LAPW 
method (with a remarkably small maximum discrepancy 
of 0.3% in structural parameters). This cross-checking 
supports the reliability of the calculations and consis- 
tency of the different methods employed. We emphasize 
that these calculations for the structure were done with- 
out the inclusion of magnetism. 

As discussed in detail in Ref. 34, there is generally a 
substantial underestimation of the pnictogen heights in 
these compounds when calculations are done in this way, 
while on the other hand in magnetic calculations LDA 
and GGA results differ, with the GGA giving much larger 
magnetic moments than experiment as well as magnetism 
that persists throughout the phase diagram in disagree- 
ment with experiment. LDA calculations done with the 
GGA structure give an intermediate state less magnetic 
than the GGA calculations but more magnetic than ex- 
periment, while in magnetic LDA calculations the As 
height is still substantially underestimated. In fact, LDA 
calculations done at the non-magnetic LDA As height 
(which agrees with the GGA As height) give the weakest 
magnetism, closest to experiment, though still overesti- 
mating the strength of the SDW. It was conjectured that 
these problems are consequences of strong spin fluctua- 
tions. 

This conjecture is supported both by comparison of 
theoretical results with experiment as well as exper- 
imental observations, such as the highly unusual in- 
creasing with T susceptibility, x(T) above the spin 
density wave ordering temperature observed in some 
compoundsi^^i^i^i Such an increasing shape suggests 
that strong magnetic correlations with the character of 
the SDW persist well above T/v; since this itinerant mag- 
netic state is driven by electrons at the Fermi energy a 
related electronic reconstruction of the Fermi surface may 
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FIG. 1: (Color online) Crystal structure of hypothetical 
LiFeSb with the relaxed structural parameters from LAPW- 
LDA total energy minimization. Note that while for clarity 
similar size spheres are used for the different atoms, from a 
crystal chemical point of view Sb"^~ anions are very large while 
Li"*" is very small. 

also be expected above T^; this may be seen perhaps in 
photoemission. This shape of x(T) persists also in the 
normal state for doped samples, where there is supercon- 
ductivity but no SDW."^^ Other evidence for strong spin 
fluctuations comes from core level spectroscopy"^— and 
transport data showing strong scattering above the or- 
dering temperature.''^ Returning to the increasing x(T), 
which continues up to high temperature, one interesting 
possibility is that this represents preformed pairs that 
can condense into either superconducting or SDW order 
and which begin forming at a very high non-observed 
temperature above which x(^) would return to a more 
normal decreasing with T shape. In any case, in this 
work where we compare the compounds, we consistently 
used the relaxed atomic coordinates from non-magnetic 
calculations. 

A requirement for a compound to be made is that the 
lattice be stable. We verified that this is the case for hy- 
pothetical LiFeSb by calculating the vibrational modes 
of the compoundii22i We find no soft or unstable modes 
and no soft elastic constants. The calculated phonon dis- 
persion curve and phonon DOS for LiFeSb are shown in 
Fig. [2l Due to larger difference in atomic weights com- 
pared to LiFeAs the phonon spectrum of LiFeSb is di- 
vided into three separated manifolds. The region of high 
frequencies (above 275 cm^^) is dominated by Li, while 
the moderate (between 200 and 275 cm^^) and low (be- 
low 150 cm"""^) frequency manifolds mainly derive from 
Fe and Sb respectively. As may be seen, all the phonon 
frequencies are safely positive and there are no optical 
phonon branches with dispersions that dip towards zero 
frequency. This shows that the PA/nmm structure of 
LiFeSb is dynamically stable. We also calculated the 
heat of formation from the elements. We obtain -0.51 eV 
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FIG. 2: Left panel: Calculated phonon dispersion curves for 
LiFeSb. Right panel: (Color online) The total and projected 
(onto atoms) phonon DOS. 

per formula unit (i.e. -49 kJ/mole formula unit), which 
indicates that the compound may be delicate but would 
at least be stable against decomposition into elements. 
Therefore we continue to discuss the magnetic and elec- 
tronic properties. 

Our main results for the electronic structure of LiFeSb, 
are given in Figs. [3l[4]and[5l which show the calculated 
band structure, electronic density of states (DOS), and 
Fermi surface, respectively. The general shape of band 
structure near the Fermi energy Ep is very similar to the 
calculated results for LiFeAs.— There are compensat- 
ing heavy hole and electron Fermi surfaces, with two elec- 
tron cylinders at the zone corner (M) and hole surfaces 
around the zone center. The hole surfaces consist of 2D 
cylindrical and small heavy 3D sections. Similar to the 
Fe-As based materials,— '^'^'^'^'^ the electron Fermi 
surface of LiFeSb may be described as two intersecting 
cylindrical sections of elliptical cross-section, with major 
axes at 90° to each other and centered at the Af point. 
We find somewhat a different hole Fermi surface struc- 
ture from LiFeAs, with only one complete hole cylinder 
at the zone center, along with two additional heavier 3D 
hole pockets. It can be seen that electron cylinders are 
more two dimensional than in LiFeAs. Also, the 2D hole 
cylinder is close in size to that of the electron cylinders, 
which may be expected to lead to nesting. Thus there 
is strong nesting of Fermi surface at the two-dimensional 
(2D) nesting vector (7r,7r). This would be expected to 
lead to an SDW state related to the M point, as in the 
Fe-As based superconductors. ^^^^i^i^i^i^ We studied 
the energetic stability of SDW state for LiFeSb directly 
using a doubled cell containing lines of Fe atoms with 
parallel spin in the Fe-Sb layers and do in fact find a sta- 
ble SDW state. Within the LDA with the LDA structural 
parameters, the local spin moment of the SDW state is 
1.12 fiB, much larger than the corresponding value (0.69 



FIG. 3: Calculated LDA band structure of LiFeSb using the 
calculated structural parameters. The Fermi energy is at 
eV. 



IIb) for LiFeAs calculated in the same way ,-22. indicating 
that LiFeSb has a more stable SDW. 

The qualitative similarity to the electronic structure 
of the Fe-As based superconductor o^^'^^i'*°''^^i"^^''*'^ is also 
evident in the DOS. The Sb p states are located mainly 
below —1.7 eV relative to Ep, and are only moderately 
hybridized with the Fe d states, indicating that the Sb 
is anionic with valence close to 3. The DOS near the 
Fermi level is dominated by Fe d states, deriving from 
the metallic Fe^"*" sublattice with direct Fe-Fe interac- 
tions, and has a characteristic pseudogap near Ep. In 
fact, Ep lies on the low energy side of the pseudogap, 
where N{Ep) is decreasing with energy but still high. 
Specifically, the value of N{Ep) = 2.2 eV"^ per Fe both 
spins is much larger than that for LiFeAs and is com- 
parable to the oxy-arsenides (e.g. N{Ep) calculated in 
the same way for LaFeAsO is 2.6 eV"^), which are the 
Fe-As compounds with higher Tc- For comparison, the 
values for LiFeAs and BaFe2As2 are 1.79 eV~^ and 1.53 
eV~^ on a per Fe basis, repectively, when calculated in 
the same wayi^ 

Within the Stoner theory, the appearance of an in- 
stability of the paramagnetic state towards itinerant 
fcrromagnetism would be determined by the criterion 
N{Ep)I > 1, where / is the Stoner parameter, with 
the typical value in Fe compounds of / 0.7 - 0.8 eV. 
Thus, the significantly higher N{Ep) in LiFeSb would 
inevitably place it closer to magnetism in general than 
LiFeAs or BaFe2As2. While the mechanism of super- 
conductivity has yet to be established, there is accumu- 
lating evidence of a connection with magnetism, and so 
chemically tuning the proximity to magnetism is a likely 
strategy for modifying the superconductivity. In gen- 
eral, the Fe-based superconductors exhibit temperature- 
induced magnetic and structural phase transitions with a 
spin density wave (SDW) character related to the Fermi 
surface nesting j^ i'^^i'^^ Superconductivity appears as the 
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FIG. 4; (Color online) Calculated total and partial electronic 
DOS for LiFeSb, on a per formula unit basis. The contri- 
bution from Li-2s state lying in deep energy range was not 
shown. The projections are onto the LAPW spheres, thus the 
Sb-5p was slightly underestimated owing to its more extended 
orbitals. 



LiFeSb 



LiFeAs 





FIG. 5: (Color online) Calculated LDA Fermi surface of 
LiFeSb in comparison with LiFeAs, shaded by band veloc- 
ity with blue as low velocity. The right panels are top views 
along the c-axis direction. 



spin density wave is suppressed by doping or pressure. 

Electronic structure calculation s^^'^^i"^°'^^''*^i'*^ show 
that all these materials have compensating small elec- 
tron and hole Fermi surfaces, with nesting between 2D 
electron sheets and heavier 2D hole sheets, which are 
separated by (tTjTt). This is associated with the SDW 
magnetic state4^^^^^^ 

Within this framework, the idea that going to heavier 
ligands may be beneficial for superconductivity is sup- 
ported by previous density functional calculations. These 



have shown that the electronic structures of the Fe-As 
superconductors are rather ionic with the exception of 
the Fe layers, which are metallic due mainly to Fe-Fe 
interactions.^^ This is different from the cuprates where 
hopping is through the O atoms in the Cu02 planes, and 
implies that the ligand (0/As) atoms play a less crucial 
role in the properties of the FeAs superconductors than in 
the cuprates. Furthermore, it has been found that there 
is a strong connection between the As position above the 
Fe plane and the magnetic properties, with higher posi- 
tions yielding stronger magnetismi^i^ This is supported 
by calculations comparing FeSe and FeTe^ and for hy- 
pothetical LaFeSbO in comparison with LaFeAsOj^ In 
both cases stronger magnetism is found in going to the 
larger ligand, which because of its size is then further 
from the Fe plane yielding narrower bands and higher 
N{Ef). The combination of a stronger SDW and higher 
N(Ep) leading to stronger spin fluctuations in general 
and in particular away from the nesting vector may be 
crucial. This is because the ordered SDW is antagonistic 
to superconductivity. In scenarios where the associated 
spin-fluctuations that couple the electron and hole Fermi 
surface sections play the main role in pairing, as dis- 
cussed in Refs. and Hi], spin fluctuations away from 
the nesting vector while not directly pairing may play a 
very important role. This is because they would compete 
with the SDW preventing long range order and leading to 
a renormalized paramagnetic state even though in mean 
held the SDW may be the predicted ground state as in the 
LDA. In any case, in these scenarios the role of doping is 
to weaken and broaden the peak in the susceptibility as- 
sociated with the nesting, destroying the SDW in favor of 
a state with spin fluctuations around the zone corner. We 
also note that in case the SDW is not destroyed by doping 
alone, it may be possible for it to be destroyed by dis- 
order yielding superconductivity in an alloy system such 
as Fei+2;Te - LiFeSb. This may be possible because the 
SDW is related to a divergence in the peak of the suscep- 
tibility, x(q), while within a spin-fluctuation mediated 
framework superconductivity will in general be related 
not to the peak value but to an integral over the Fermi 
surface, i.e. a high average value over some region of 
the zone, specifically the region for which the wavevector 
connects the electron and hole Fermi surfaces^^*^ Also 
near divergences will be pair breaking for superconduc- 
tivity. This may also partly explain why superconduc- 
tivity in these phases is relatively robust against alloying 
with Zn or Co in the Fe planes^^i^i^ even though in 
an unconventional superconductor scattering, including 
non-magnetic scattering, is pair breaking. Thus, even if 
the magnetic ground state cannot be destroyed by doping 
in Fei+2;Te or LiFeSb (supposing that this can be synthe- 
sized) it may be destroyed in favor of superconductivity 
in the solid solution between these two compounds. 

In any case, our results show that, if it can be syn- 
thesized, LiFeSb will have electronic and magnetic prop- 
erties closely related to those of the Fe-As based super- 
conductors, and in particular will show a qualitatively 



5 



similar Fermi siirfacc structure and tendency towards an 
SDW state. In comparison with LiFeAs, it will have 
a higher N{Ef) and a stronger SDW. This may favor 
higher critical temperatures, at least within a scenario 
with interband pairing mediated by spin fluctuations as- 
sociated with the Fermi surface nesting. In addition, this 
material is found to be dynamically stable, evidenced by 
the absence of any unstable phonon modes. It is also 
worth noting that this solid solution contains no elements 
as toxic as As. As such it would be of considerable in- 
terest to attempt synthesis of this compound or its solid 



solution with Fei+a;Se, LiFeAs and especially Fei+^Te. 
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